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1- A Brief History of Research on the 
Aeroacoustics of Rotating Blade Machinery at 
Langley 

Langley has been in the forefront of research on aircraft noise generation and 
prediction since the 1930's [SI]. Some of the early researchers in this subject at 
Langley were Deming, Garrick, Watkins, Hubbard and Maglieri. This research was 
aided greatly with the introduction of the high speed digital computers in the 1960's 
that were used both in acoustic experiments as well as in the development of 
noise prediction codes. About the same time, the pace of research in the field of 
aeroacoustics in universities, industry and government agencies accelerated due 
to the public objection to aircraft noise both in the cabin and around the airports. 
Modern digital computers have made it possible to use advanced mathematical 
models to study the acoustics of propulsion systems with realistic geometry and 
operating conditions. NASA scientists and engineers have contributed significantly 
to aircraft noise reduction both by theory development and the use of experimental 
methods. 

In what follows, we will discuss the history of theoretical model development at 
Langley with which the author was involved during the period 1975 to present. We 
will confine ourselves to subjects of relevance to aircraft with open rotors as their 
propulsion system. This would also involve the calculation of installation effects 
and the noisefoot prints at various flight conditions. 

The introduction of the concept of the acoustic analogy (AA) by Lighthill in 1952 
[S2] was one of the most important advances in aeroacoustics which has resulted 
in current high level of aircraft noise prediction technology. In 1969, Ffowcs 
Williams and Hawkings presented an extension of Lighthill's jet noise equation to 
the noise from turbulence in the presence of moving surfaces appropriate for 
studying the acoustics of rotating machinery - helicopter rotors, propellers, fans 
and open rotors [S3]. This extension is known as the Ffowcs Williams-Hawkings 
(FW-H) equation. 

Let a moving surface be described implicitly by f{ x, t) = 0 such that vf = n, n is the 
unit outward normal. This assumption implies that f > 0 outside the moving surface 

(see Figure 1). 
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f(x,t)< 0 


Figure 1- The definition of the moving surface implicitly as f(x, t ) = 0. Note 
that v f - n where n is the unit outward normal to the surface. 


The Ffowcs Williams- Hawkings (FW-H) equation is 

(3 (3 ( 3 ^ 

□ 2 P' = - [po V„ <S(/)] - — [p n, <?(/)] + — — mf) Tij ] 

Ot OXi OXjOXj 


( 1 ) 


Here and elsewhere in the paper, we will use the summation convention on the 
repeated index. In this equation n 2 is the wave or D’Alembertian operator in three 
dimensional space. Also p' = c 2 p' = c 2 (p-p 0 ), c and p 0 are the speed of sound 
and density in the undisturbed medium, respectively. Note that p' can only be 
interpreted as the acoustic pressure if p'/p 0 « J \. The symbols v n , p and 

Tjj - p Uj Uj-crij + (p' - c 2 p')6ij are the local normal velocity of the surface, the 
local gage pressure on the surface (in fact p-p 0 ), and the Lighthill stress tensor, 
respectively. In the definition of the Lighthill stress tensor, cr,y is the viscous stress 
tensor and <5,y is the Kronecker delta. The Heaviside and the Dirac delta functions 

are denoted H(f) and 6(f), respectively. In the second term on the right of eq. (1), 
we have neglected the viscous shear force over the blade surface acting on the 
fluid. 
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The three inhomogeneous source terms on the right of Eq. (1) are known as the 
thickness, loading and quadrupole sources, respectively. Note that the thick- 
ness and loading sources are surface sources while the quadrupole source is 
regarded as a volume source although it can be analytically broken into volume 
and surface sources [21, 23, 27, 31, 36]. This latter point of view can explain 
several (old and new) mechanisms of noise generation in rotating blades some of 
which have not been studied in detail, e.g., flapping wakes in the vicinity of trailing 
edge. 

We note that we have artificially converted a nonlinear problem of noise generation 
by a moving surface to a linear problem by using the acoustic analogy. All the 
nonlinearities are lumped into the Lighthill stress tensor which is assumed known 
from near field aerodynamic calculations. When we started working on helicopter 
and propeller noise in the early seventies, because of the limitations of digital com- 
puters, the most we could expect from aerodynamic calculations was the blade 
surface pressure. For this reason, using some physical reasoning, we neglected 
the quadrupole volume sources in FW-H equation and concentrated on develop- 
ment of formulations for the prediction of thickness and loading noise. Later on, as 
computers became more powerful, we included quadrupoles in our noise predic- 
tion. There was, however, another theoretical advance which led to the use of 
purely surface sources to which we will turn next. 

It was Ffowcs Williams himself who proposed to use a penetrable (porous or perme 
able) data surface to account for nonlinearities in the vicinity of a moving surface. 
We again assume that the penetrable surface defined by f(x, t) = 0 and the fluid 
velocity is denoted by u. The FW-H equation for penetrable (permeable, 
porous) data surface, FW-H pds , is: 



(2) 


We have used the following notations in the above equation: 



(3) 


Lj — p 5 jj Hj + p Uj {u n v n ) 


(4) 
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where 8,j is the Kronecker delta. As in the case of eq. (1), in the first term on the 
right of eq. (4), we have neglected the viscous shear force over the data surface 
acting on the fluid exterior to the surface. The philosophy behind using FW-H pds is 

to locate the data surface f = 0 to enclose a moving surface, in such a way that all 
quadrupoles producing non-negligible noise are included within this surface. There- 
fore, no volume integration of the quadrupoles outside the data surface is neces- 
sary. High resolution CFD calculation is performed in the near field region 
(including turbulence simulation if broadband noise prediction is required). The 
data surface used for the acoustic calculation should be located within the region 
of high resolution CFD computation. As discussed in Section 3, the optimal loca- 
tion of the data surface, particularly when vortices cross the surface, is still the 
subject of research. One would like this surface to be as small as possible, 
because of the computer intensive nature of aerodynamic and turbulence simula- 
tion that require fine grid sizes and small time steps. 

We have developed several equivalent solutions of the FW-H equation in the time 
domain (called formulations in this report) at NASA Langley that are used to 
develop two efficient codes for the prediction of the noise of rotating machinery. In 
the next section we present a discussion of these formulations, computer codes 
based on them and their capabilities. Note that the Fast Scattering Code (FSC) 
does not depend on these formulations but is included in the next section because 
we will find it useful for open rotor noise prediction when scattering of the open 
rotor generated noise by fuselage and wing surfaces becomes important. 
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2- Noise Prediction Codes for Open Rotors 
Developed at NASA Langley and Their 
Capabilities 

The author has derived several solutions of the FW-H equation that are useful for 
efficient numerical computation of the noise of rotating blades. We mention here 
three formulations in the time domain valid in the near and far fields that 
include only the thickness and loading sources of the FW-H equation. In chronologi- 
cal order, these are: 

1- Formulation 1 [1,8, 52]: This was developed for surfaces moving at subsonic 
motion if used with the Doppler singularity. It is valid for surfaces moving at tran- 
sonic and supersonic speeds if at no time during the emission process the condi- 
tion M r - 1 occurs and the collapsing sphere form of the formulation is utilized [1]. 
Here M r is the Mach number of the source in the radiation direction. Note that in 
this section, the Mach number is always based on the speed of sound in the undis- 
turbed medium. In the case of transonic and supersonic speeds, the multiple emis- 
sion times of the source must be accounted for in the acoustic calculation. Formula- 
tion 1 is the simplest of our three formulations distinguished by an observer time 
derivative before one of its two surface integrals. This time derivative must be 
evaluated numerically which can be a source of numerical error. 

2- Formulation 1A [7-9, S4, 52]: This result was obtained by taking the observer 
time derivative of Formulation 1 inside the integral analytically and, thus, eliminat- 
ing the need for taking the derivative numerically. Although, in principle, this result 
can be obtained by elementary calculus (not recommended), the simplest method 
of deriving it is using an identity of generalized functions for the free space Green's 
function of the wave equation at the beginning of the derivation [8]. An immediate 
consequence of eliminating the observer time derivative of Formulation 1 is the 
increase in the complexity of the formula and the increase in the order of the 

Doppler singularity from (1 -M r ) _1 to (1 - M r )~ 3 . The increase in the complexity of 
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the formulation does not cause any concern for numerical implementation on the 
computer. But the increase in the order of the singularity increases the sensitivity 
of the result near M r = 1. However, in practice, we must switch to Formulation 3 in 
the acoustic code when this situation occurs. The conditions for the validity of 
Formulation 1A is similar to Formulation 1 and is used in the forms with Doppler 
singularity (subsonic sources) as well as collapsing sphere (transonic and super- 
sonic sources). An important use of Formulation 1 A is in the acoustic code based 
on the FW - H pds equation. 

Formulation 3 [13, 16, 17, 20]: This complicated result is very difficult to derive 
and requires advanced knowledge of differential geometry, general tensor analysis 
and generalized functions. The formulation is valid at all speed ranges of the 
source surface - subsonic, transonic and supersonic. However, because of its 
complexity, it is recommended for use it of the parts of the surface for which the 
condition M r = 1 occurs. For an open rotor with blades that have sharp leading and 
trailing edges in the region that run supersonically, the blade sweep can be 
specified in such a manner that no singularities would appear in the acoustic 
pressure signature when this formulation is used. Otherwise, one can expect the 
integrable logarithmic singularity in the acoustic pressure [22, 27]. However, it 
should be mentioned that the appearance of this singularity is mathematical and 
can be removed by the inclusion of the quadrupole term of the FW-H equation as 
shown by Farassat and Myers [38]. 

An explanation of why Formulation 3 is so complicated is given in [46]. An impor- 
tant test of this formulation is to compare the predicted acoustic pressure signa- 
tures and spectra from a subsonic rotating blade exclusively from Formulations 1A 
and 3 under identical input data. This computation has been performed at NASA 
Langley, documented but not published. The agreement of the predicted results 
from the two formulations was excellent for both the thickness and loading noise. 
Furthermore, this agreement also verified that both these formulations were coded 
correctly in the Langley developed code ASSPIN (Advanced Subsonic - Super- 
sonic Propeller Induced Noise). 

Some old and obsolete NASA Langley acoustic codes used Formulation 1 for high 
speed propellers [6]. We currently use only Formulations 1A and 3 in our noise 
prediction codes such as ASSPIN discussed in the next section. The best article to 
read about how to implement these formulations most efficiently to predict the 
noise of rotating blade machinery is [20]. 
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We now turn to the description of two computer codes for open rotor noise predic- 
tion. Note that our codes are valid for an observer in the near and far fields. 
The third code we discuss below is the Fast Scattering Code (FSC) which is recom- 
mended for the prediction of the effect of the scattering of the open rotor noise by 
fuselage and wings of an aircraft. This phenomenon includes the important shield- 
ing effect of noise by the rigid surfaces of the aircraft known also as the installation 
effect. 

■ 2.1- Advanced Subsonic-Supersonic Propeller Induced Noise 
(ASSPIN) Code 

This code was developed in mid-nineteen eighties originally for a single row of 
open rotors known as propfans that operated at supersonic tip speed at cruise 
condition [13, 16, 17, 20]. It can also calculate the noise of a contra-rotating pro- 
peller with equal number of blades in each row. However, the present code can 
form the core of a much more capable open rotor noise prediction code by the 
addition of pre- and post-processors. There is an excellent Users Manual written 
by Dunn and Tarkenton [S5]. This code uses Formulation 1A and 3 in such a way 
that the latter formulation is used for panels on the blades with near sonic speed. 
Farassat, Padula and Dunn [20] have described the structure and the algorithms of 
ASSPIN in great detail with some comparison of the predicted and measured 
acoustic data for propfans. Other publications of Langley with more examples of 
application of the code are [30, 34, 35]. One of the most significant tests of 
ASSPIN was against the frequency domain code of Donald B. Hanson with identi- 
cal blade geometry, operating conditions and blade aerodynamics. The results 
with good agreement were published in the same issue of the AIAA Journal in 
1992 [35, S6] 

The ASSPIN code was written in FORTRAN 77 and developed originally on a 
CDC-7600 mainframe computer of NASA Langley. It runs very efficiently on a 
modern scientific desktop computer. Noise footprint computations can be per- 
formed routinely in a matter of few hours and even the animation of the footprint 
during the takeoff condition can be simulated in a reasonable time. Other than the 
fact that the nonlinear effects are not included in this code, it is one of the most 
robust, accurate and verified codes in the time domain in U.S. and, to our best 
knowledge, in the world. 

□ 2.1 .1 - Input Data to ASSPIN 

1- Hot blade geometry- several options are available to users 

2- Number of blades in each row 

3- Inflow angle to the rotor axis 

4- Unsteady blade surface pressure 
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2.1 .2- Output Data From ASSPIN 

1- Acoustic pressure time history for a stationery observer with respect to the 
ground or moving with the aircraft at uniform speed 

2- Acoustic pressure spectrum for an observer moving with the aircraft at uniform 
speed or the acoustic spectrum for contiguous segments of time for a stationary 
observer (nonperiodic time history). 

■ 2.2- Ffowcs Williams-Hawkings with Penetrable Data Surface 
(FW - H pds ) Code 

It has been known since the early days of high speed propeller noise research that 
nonlinear propagation effects in the near field are important. To include the nonlin- 
earities, as discussed in Section 1, the use of a prediction code based on the 
FW-Hp ds is an obvious choice. Since all the commercial airliners will fly at sub- 
sonic cruise speed in the foreseeable future, Formulation 1A is the appropriate 
solution of the FW-H pds equation to use. For an aircraft in steady rectilinear flight 
(the most common situation of interest), there is a simple closed-form formula 
based on the Garrick triangle construction for the computation of the retarded time 
in the acoustic calculations. Thus, the development of a noise prediction code 
based on FW-Flp ds equation using Formulation 1A for this case is a relatively 
simple task. Nark and Farassat have developed such a code at Langley and have 
used it for ducted fan noise prediction [50]. It can be used for open rotors with 
some additional pre- and post-processors. 

When turbulence and wakes from rotors cross the data surface, spurious acoustic 
signals appear in the acoustic signature which reduce the accuracy of this method. 
We will discuss in Section 3 our plans to alleviate this shortcoming. 

□ 2.2.1- Input Data to FW -H pds Code 

High resolution unsteady CFD data (static pressure, fluid density and velocity) on 
the data surface must be specified. 

□ 2.2.2- Output Data From FW - H pds Code 

1- Acoustic pressure time history for a observer stationary with respect to the 
ground or moving with the aircraft at uniform speed 

2- Acoustic pressure spectrum for an observer moving with the aircraft at uniform 
speed or the acoustic spectrum for contiguous segments of time for a stationary 
observer (nonperiodic time history). 
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■ 2.3- The Fast Scattering Code (FSC) 

This code was developed at NASA Langley by Dunn and Tinetti based on the 
equivalent source method (ESM) [S7, S8]. The theory behind the FSC is 
explained well in these references. It is based on the equations of time harmonic, 
linearized acoustics and the ESM is used for solving an exterior Helmholtz equa- 
tion as a boundary value (BV) problem for a surface moving at uniform rectilinear 
speed. The incident sound source, e.g., the open rotor noise, is Fourier analyzed 
and the FSC is then used to compute the scattered noise for each frequency com- 
ponent. The ESM is a much more efficient method of finding the scattered noise 
than other well known techniques such as the boundary element method (BEM). 
The advantage comes from the fact that in ESM one does not have to deal with 
singular integrals (as in the case of BEM) which require elaborate and time consum 
ing algorithms for their evaluation. For this reason, the use of ESM in the scatter- 
ing problem results in a simple and efficient code. 

The FSC has been used to study the shielding effect of the engine noise by a 
wedge shaped airframe for which some experimental acoustic data were available 
[S9]. It was found that the agreement between the measured and predicted acous- 
tic data was good. It is clear from this example that the shielding effect of fuselage 
and wings can be an important method of aircraft engine noise control. This effect 
can only be studied and quantified by utilizing the FSC. Note that the FSC pro- 
vides some simple models of incident engine noise. 

The ESM requires the inversion of a large matrix whose size and memory require- 
ment increase with frequency co. Typically, the memory needed for the code is 

proportional to (jj 4 and the computational time is proportional to The proportion- 
ality constants are functions of the scattering surface area. It can be shown that on 
today's scientific desktop computers with one processor, either the size of the 
scattering surface must be kept small or the incident frequency must be assumed 
to be low [S8]. Fast multipole method (FMM) can be used to increase the upper 
limit of the frequency for the same scattering surface by a factor of 3 to 4 [SI 0]. 
However, even with the use of FMM, one is not able to calculate the scattering of 
noise from a large aircraft, say the size of Boeing 777, and an incident frequency 
of 10 khz on a scientific desktop computer. The FSC has been parellelized to run 
on multiprocessor machines to achieve this goal. 

The FSC is written in FORTRAN and can be operated in batch mode. Version 2.0 
is written for a single processor. It can also be driven by a graphical user interface 
(GUI) written in C++ based on WXWidgets. Two parallel versions of the FSC have 
been developed at Langley - v3.1 for shared memory machines, and v3.2 for 
distributed memory machines. 
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The FSC was developed for the UNIX and LINUX operating systems but it can 
also be adapted for Windows and Macintosh environments. A Users Manual is 
available for various versions of the code [S8]. 

We believe that FSC is the most advanced and efficient aeroacoustic code avail- 
able today for computation of the scattering of incident sound by aircraft surfaces. 
With parallelization to run on multiprocessor machines, it can fully meet the needs 
of the aircraft and engine industry for product design, evaluation and noise reduc- 
tion. 

□ 2.3.1- Input Data to FSC 

1- Excitation frequency and freestream thermodynamic variables 

2- Numerical parameters for determining the appropriate grid generation and solu- 
tion strategies 

3- The scattering surface geometries (wings, fuselages, nacelles) and liner 
admittances 

4- Local flow variables (density, speed of sound, and Mach number vector). Code 
options include: a- no flow, b- uniform flow, c- FSC generated small perturbation 
compressible flow, and d- user supplied flow 

5- Complex values of incident acoustic pressure and acoustic velocity at FSC 
requested locations in space and time. Code options for incident sound include: a- 
simple point monopole(s) and dipole(s), b- FSC generated engine noise from 
nacelle alone run, and c- results from external noise codes capable of producing 
FSC type inputs. 

6- Observer position(s): Specify arbitrary observer points or request field calcula- 
tions on spheres, cylinders, rectangular volumes surrounding the scattering sur- 
face. An option for requesting take-off or flyover footprints is provided. 

□ 2.3.1- Output Data From FSC 

Complex values of the instantaneous acoustic pressure, velocity and intensity at 
the desired observer positions or the noise footprint regions. 
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3- Research and Development Plans 

It can be seen from the above presentation that Langley researchers have devel- 
oped the aeroacoustic technology to study the open rotor noise generation and 
reduction. The available codes are able to meet much of the need of the aircraft 
and engine industry for product design, evaluation and noise control. We feel that 
there are areas that are in need of further improvements. We also feel that the 
aircraft and engine industry have specific problems of interest that may require 
further extensions of the capabilities of our codes. We are interested to learn of 
these problems. We present below some of our research and development plans 
which will increase the capabilities and usefulness of our codes. 

■ 3.1- Proposed Plans for ASSPIN Code 

ASSPIN is in an advanced stage of development and we suggest that 
this code be used for studying free open rotor noise. The reason is that 
to use the FSC with ASSPIN output as the incident sound source, we 
must be able to calculate the gradient of the acoustic pressure on all 
the scattering surface. The most numerically efficient method of finding 
the gradient is by an analytic method. This can be done rather easily for 
Formulations 1 and 1A [54] but it is next to impossible for Formulation 3 
because of its analytic complexity and the number of terms. However, 
since the ASSPIN code was developed in nineteen eighties, there is a 
need for upgrading some of its capabilities. We recommend the follow- 
ing improvements to the ASSPIN code: 

1- Improvements of the blade geometry specification- The blade geometry of 
the modern open rotors is a complex three dimensional curved surface. We 
require the calculation of various curvatures of the surface to use in Formulation 3. 
This means that we need a surface specification with smooth second derivatives. 
This task was relatively easy when NACA 2D airfoils were used by some open 
rotor blade designers because analytic expressions for airfoil coordinates were 
available. However, today's blade designers use advanced airfoils optimized by 
CFD codes for their performance. We need to use sophisticated three dimensional 
blade surface specification, e.g., NURBS, which will give us surface curvature 
information for noise prediction. We have made some progress in this direction but 
more work needs to be done in this area. This feature will be added as a pre-pro- 
cessor to ASSPIN. 
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2- Upgrading ASSPIN to compute the noise of contra-rotating open rotors 
with unequal blade numbers and geometry in each row- Currently, ASSPIN 
code can handle contra-rotating open rotors with equal number of blades on each 
row and identical geometry. To treat the general case of unequal number of blades 
on each row, the acoustic pressure signatures of the interaction tones must be 
computed. This involves changes to the observer time step and period specifica- 
tions in ASSPIN code itself. Furthermore, the part of the code for the computation 
of the acoustic spectrum must be modified to account for the interaction tones. 
This feature will also give us the ability of predicting the noise of a contra-rotating 
open rotor with two rows of rotors with very different geometries (blade numbers 
and blade designs). 

■ 3.2- Proposed Plans for FW - H pds code 

We have not used this code for open rotor noise prediction and what we are propos 
ing is all new and requires further research. However, for two reasons we feel that 
work on this methodology must proceed with urgency because of the important 
capabilities it will provide us. First, we can include the near field nonlinear propaga- 
tion effects which can be important in noise prediction. Second, this is the only 
open rotor prediction method that can be linked easily to the FSC to study the 
shielding effect and noise footprint changes by scattering from fuselage and wings 
(see Subsection 3.3, Item 3). We, therefore, propose the following research plans: 
1- Removal of spurious signals from the computed acoustic pressure signa- 
ture when eddies and wakes cross the data surface - When turbulent eddies or 
the wakes from open rotors cross the data surface of the FW - H pds code, spurious 
signals appear in the acoustic pressure signature. This is a major cause of con- 
cern for the accuracy of this method. We now outline here a method to alleviate 
the spurious signals. Our discussion refers to Figure 3 below. 
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Figure 3- Sketch of the data surface for FW - H pds code and the use of equiva- 
lent sources to generate the null field inside the data surface at control 
points. Here we show one equivalent source and one control point for 
simplicity. 
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We first note that in the application of this code, we specifically ignore the 
quadrupoles outside the data surface because specifying these sources is not 
practical. We only use surface sources in this method. When the rotor wakes cross 
the data surface, we have a fluctuating quadrupole field in the exterior of the data 
surface which generate the spurious signals. We have one important piece of 
information about the quadrupole field from the surface data. The sum of the contri- 
butions of the quadrupoles and the surface sources inside the data surface is identi- 
cally zero (null field) for all time. Therefore, we can use some discrete sources 
(generally monopoles that we call equivalent sources) with unknown strengths in 
specified exterior or interior locations near the area of wake crossing to produce 
null field at specified control points inside the data surface. We must have, at a 
minimum, as many control points as the number of equivalent sources to be able 
to find the unknown strengths of the equivalent sources from a system of linear 
equations. Once the strengths of equivalent sources are found, the noise at any 
observer position in the exterior of the data surface is found from the surface 
sources as well as the equivalent sources and the spurious signals should disap- 
pear. 

The problem of finding the source strengths is almost certainly ill-posed of the 
common variety seen, for instance, in solving the scattering problem using equiva- 
lent sources (as in FSC). The handling of such ill-posed problems is understood 
quite well and we do not expect undue difficulties to find the strengths of the equiva- 
lent sources. 

2- Development of an algorithm to transfer CFD data to the data surface for 
acoustic calculations - This is an important recurring problem in aeroacoustics. 
We ask for high resolution CFD data for noise calculations and invariably the CFD 
grids do not coincide with the grids on the data surface for acoustic calculations. 
Using a surface generated by the CFD grids for acoustic calculations generally 
does not work because we need some surface geometry information (outward 
normal, etc.) that cannot be obtained easily from such a surface. The accurate 
extrapolation of high resolution CFD data to the data surface presents a challeng- 
ing problem that needs to be solved to make the FW - H pds code work satisfacto- 
rily. Such an algorithm must be implemented as a preprocessor for the code. 
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■ 3.3- Proposed Plans for the Fast Scattering Code (FSC) 

The FSC is in the process of intensive development for NASA Langley by Dunn 
and Tinetti. It has not been used for calculation of scattering of sound for an open 
rotor as the incident source of sound. However, we feel that the scattering of open 
rotor noise from aircraft fuselage and wings must be studied in the design stage of 
the aircraft for noise control. We propose the following developments and research 
plans for the FSC: 

1- Search for efficient and accurate method of storage and inversion of large 
and dense matrices on multiprocessor machines - For multiprocessor 
machines, the selection of an efficient storage and the inversion method of large 
and dense matrices that ESM requires is a challenging problem. The accuracy of 
the results and the execution time are both dependent on the method of storage 
and inversion of the matrix. One criterion of the search for an algorithm is reaching 
a better balance between memory and time for the very large cases. This work will 
be performed by Dunn and Tinetti for NASA Langley. 

2- Development of a method for finding the pressure gradient of the incident 
sound at all points of scattering surface - To satisfy the boundary condition on 
the scattering surface (equality of the normal velocities of incident and scattered 
acoustic pressures), one needs to know the gradient of the incident sound. This is 
a very difficult problem. In the case of using ASSPIN for open rotor noise predic- 
tion, we have not been able to solve this problem yet. But for the case of FW - H pds 
code used to simulate the incident sound for open rotors, we have been able to 
derive two analytic results for the acoustic pressure gradient based on Formula- 
tions 1 and 1A [54]. These are known as Formulations G1 and G1 A with similar 
characteristics as Formulations 1 and 1A. Either result should be coded as a post- 
processor of the FW - H pds code and thus provide the input data to the FSC. 
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